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Abstract 

A  high-pressure  chamber  is  used  to  investigate  and 
further  enhance  our  knowledge  and  physical 
understanding  on  effects  of  thermodynamical  subcritical- 
to-supercritical  transition  of  ambient  condition  ^on 
cryogenic  liquid  injection  using  two-dimensionaTjJlaman 
scattering.  Pure  liquid  N2  is  injected  into  N2.  The 
injector  is  a  50$fmicron  diameter  straight  hole  having  a 
long  length-to-diameter  ratio  of  100.  The  optical  setup 
uses  a  pulsed  Nd:Yag  laser  frequency-doubled  to  532 
nm.  Difficulties  arise  with  optical  breakdown  of  the  N2 
molecules  in  drops  and  ligaments  by  local  focusing  of  the 
laser  beam  dominating  the  Raman  signaf;  particularly  at 
sub-  and  near-critical  regions.  The  severity  of  this 
problem  is  reduced  by  stretching  the  laser  pulse  width 
using  a  double-loop  design  with  mirrors  and  beam 
splitters.  Careful  and  painstaking  alignment  is  needed  to 
take  advantage  of  this  pulse-stretcher  design.  Two- 
dimensional  images  are  taken  near  the  injector  and  results 
interpreted  in  terms  of  density  plots.  At  subcritical 
ambient  conditions  a  small  number  of  images  are  needed 
for  averaging  and  strong  Raman  signal  is  obtained. 
Substantial  reduction  of  this  signal  is  observed  at 
supercritical  condition  due  to  jet  radial  growth  and 
mixing.  Radial  intensity-surplus  profiles  are  broadened 
by  the  interaction  of  the  laser  and  the  jet  interface.  The 
corrected  profiles  show  a  tendency  towards  the  similarity 
profile  seen  for  the  gaseous  turbulent  jets  near  and  above 
the  critical  point.  Initial  growth  rate  of  the  jet  as  judged  by 
the  Raman  signature  shows  tendency  towards  gaseous 
turbulent  jet  when  chamber  pressure  is  raised  from  near 
the  critical  value  and  beyond.  However,  the  difference 
remains  quite  distinct  being  nearly  half  the  growth  rate 
reported  for  gaseous  jets.  The  initial  growth  rate  is  in 
reasonably-good  agreement  with  our  earlier 
measurements  using  shadowgraphy  if  twice  the  FWHM 
of  the  normalized  intensity  plots  are  used.  Results  are 
compared  with  gaseous  jets  and  conclusions  are  drawn 
regarding  effects  of  supercritical  condition  on  the  jet 
structure. 

Introduction 

It  has  been  recognized  that  there  were  advantages  in 
operating  most  power  and  thrust  producing  engines  at 
higher  combustion  chamber  pressures  and  temperatures 


for  enhanced  efficiency  and/or  output  thrust  levels.  In 
this  paper,  however,  the  main  interest  is  on  the  cryogenic 
liquid-propulsion  rockets.  Progress  in  material  science 
and  technological  achievements  in  designing  vessels  to 
withstand  high  pressure  and  temperature  conditions  bring 
us  close  to  harness  the  theoretical  gains  in  high  pressure 
and  temperature  engine  operation.  However,  there  are 
some  new  issues  that  arise  as  a  result  of  the  elevated 
chamber  conditions  which  must  be  understood  for  better 
engine  design  and  performance  predictability.  One  such 
issue  is  the  realization  that  in  cryogenic  rockets, 

particularly  under  full  throttle  or  take-off  conditions, 
liquid  propellants  may  find  injected  into  a  \ 

thermodynamically  supercritical  condition.  Our 
understanding  of  the  injection  under  these  conditions  has 
been  limited  until  recently  when  progress  was  made  in 
both  experimental  and  theoretical/computational  areas, 
see,  for  example,  Mayer  et  al  [1,  2],  Chen  and  Sui  [3], 

Woodward  and  Talley  [4],  Harstad  and  Bellan  [5], 

Delplanque  and  Sirignano  [6],  Oefelein  and  Yang  [ 7 7, 
and  Chehroudi  et  al  [ 8 ,  9]. 

In  most  experimental  works,  initial  attempts  were  to 
visualize  the  drop/jet  formation  and  evolution  using  high 
speed  photography  techniques.  Valuable  information  was 
gathered  and  analyzed  using  such  methods  and  reported  in 
the  aforementioned  publications.  In  an  attempt  to  acquire 
a  more  detailed  and  quantitative  information  from  fluid 
jets  injected  into  supercritical  condition,  several  laser 
diagnostics  were  evaluated  and  one  of  the  most  suitable 
candidate  was  identified  to  be  the  spontaneous  Raman 
scattering,  see  Anderson  et  al.  [10].  Raman  scattering  was  ^ .  fYAj£ 
used  by  Woodward  and  Talley  [4]QClJex]Jcer  et  ~ alflljr  £* 

and  Oschwald  and  Schik  [12],  Oschwald  et  al.  [13]  to  1 
investigate  the  structure  of  the  nitrogen  fluid  jet  under 
supercritical  condition.  One  of  the  fundamental 
difficulties  in  extracting  true  density  mformation  from 
the  Raman  effect  was  the  uncertainty  in  the  functional 
dependence  of  the  Raman  scattering  cross  section  on 
temperature  under  these  conditions.  Raman  cross 
section  in  liquids  is  larger  (by  1.73  for  LN2  with  n=  1.2) 
than  the  gaseous  phase  by  a  correction  factor  that 
depends  on  the  refractive  index  (n).  Not  being  strictly  a 
liquid,  an  estimate  of  10%  was  given  by  Decker  et  air 
[11]  under  supercritical  condition  for  N2  at  6  MPa.  The 
other  practical  problem,  enhanced  by  the  lens  effect  of  the 
liquid-like  jet,  is  the  onset  of  stimulated  Raman  scattering 
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(same  wavelength  as  the  spontaneous  Raman)  and 
plasma  formation  (broad-band  radiation)  interfering  with 
the  linear  (power-to-density)  spontaneous  Raman. 

-  ^nderson  et  al  [10]  attempted  to  reduce  this  problem  by 
a  design- of  a  passive  pulse  stretcher  in  the  path  of  the 
incident  beam  (measured  O2  at  580nm  with  Nd:YAG 
frequency  doubled  at  53^Wi,  30C^nJ;  from  J^s  to^2}ns 
FWHM),  and  Decker  et'al  /7  7  through  ‘  usage"  of  a 
double-oscillator  Excimer  laser  producing  two  back-to- 
back  lowered-power  2-D  pulses  (measured  N2  at  382(im 
with  XeF  laser  at  35  film,  15GmJ;  total  pulse  widui  of 
4Qhs).  Oschwald  and  Schik  [12]  employed  a  l.5(W 
continuous  wave  Ar+  laser  at  48&{un  and  detected  Raman 
signal  in  a  N2  jet  at  550.6  nm  scattered  from  within  the 
beam  focus  point  under  supercritical  condition.  No  special  « ^ 
provision  was  used  for  removal  of  the  stimulated  Raman,  tjtfP 
In  these  experiments,  Decker  et  at [11]  observed  a^f r 
smooth  transition  from  the  high-density  core  of  an  N2- 
into-N2  jet  towards  the  gaseous  outer  region;  that  is,  there 
was  no  distinct  phase  interface.  The  density  radial  profile 
was  never  a  top-hat  even  as  close  as  a  distance-to- 
diameter  ratio  of  2.5.  Oschwald  et  al  [13]  performed 
injection  of  N2  into  N2  under  supercritical  chamber 
pressures  (Pr=1.17  and  1.76)  and  at  three  different  initial 
injection  temperatures  (100,  1 18,  and  140  K).  Although 
scatter  can  be  seen  in  their  data,  the  centerline  density 
axial  profile  at  Pr=1.17  relaxed  to  the  chamber  density 
faster  the  higher  the  injection  initial  temperature.  At  the 
Pr  of  1.76,  however,  insensitivity  of  the  centerline 
density  to  the  initial  temperature  was  observed.  The 
corresponding  centerline  temperature  axial  profile  at  Pr  of 
1.17,  calculated  through  an  equation  of  state,  showed  a 
constant-value  region  slightly  away  from  the  injector, 
being  close  to  the  temperature  at  which  the  specific  heat 
as  a  function  of  temperature  reaches  maximum  (called  it 
Tjfiax).  This  was  observed  only  for  the  cases  when  the 

was  less  than  the  Tmax; 


initial  injection  temperature 


-  Under  this  condition,  the  centerline  temperature  increased 

downstream  until  it  reached  T|na^.  Therefore,  most  of  the 
heat  transferred  was  used  to  expand  the  jet  (reduce  the 
density)  and  not  to  increase  the  temperature.  This  general 
behavior  was  not  maintained  with  the  same  degree  of 
contrast  as  chamber  pressure  was  raised  higher  (at  Pr  = 
1.76)  away  from  the  critical  value  of  N2.  Based  on 
limited  tests,  they  also  concluded  that  the  thermodynamic 
state  of  an  injected  supercritical  fluid  is  of  prime 
importance  in  the  jet  growth  and  not  the  injection 
velocity  or  momentum. 


Experimental  Setup 

High  Pressure  Chamber  Facility 


Figure  1.  The  liquid  N2  from  a  dewar  is  used  to  cool 
and/or  liquefy  the  injectant  passing  through  the  cryogenic 
cooler  prior  to  injection.  In  addition  to  N2,  a  branch  is 
connected  for  possibility  of  introducing  other  gases  or 
mixtures  of  gases  as  a  chamber  medium.  The  mass  flow 
rate  of  the  injectant  is  regulated  and  also  measured  via  a 
mass  flowmeter,  and  a  precision  micrometer  valve.  For 
more  details  the  reader  is  invited  to  refer  to  Chehroudi  et 
al.  [8,  9].  The  injection  in  this  study  is  through  a  sharp- 
edged  50  mm  long  stainless  steel  tube  with  a  1.59  mm 
(1/16”)  diameter  and  a  508  micron  (0.0203;mtier  hole 
io  of  100).  With  the  Rey 


(length-to-diameter  ratio  of  100). 
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The  stainless  steel  chamber  in  Fig*’  1  can  withstand 
pressures  and  temperatures  of  up  to  137  atm  (2000  psi) 
and  473  K,  respectively,  and  lias  two  facing  circular 
sapphire  windows  and  two  UV-grade  side-mounted  slot¬ 
shaped  quartz  windows  for  optical  diagnostics  and  laser 
illumination  of  the  chamber.  Transport  of  the  liquid 
through  the  lines  and  the  injector  is  achieved  as  shown  in 


Figure  1.  Schematic  diagram  of  experimental  setup  for  sub-  to  supercritical  jet  Injection 

number  ranging  from  25,000  to  70,000,  the  entrance 
length  of  50  to.  100  is  needed  for  fully-developed 
turbulent  flow^  (§ee  Schlichting  [ 74j£rThe  length  is 
therefore  long  enough  to  ensure  fully-developed 
turbulent  pipe  flow  at  the  exit  The  rig  is  fully 
instrumented  with  thermocouples,  pressure  gages, 
transducers,  and  mass  flowmeters  at  the  locations 
indicated  in  Fig.  1. 

Optical  Setup 

An  Nd:Yag  pulsed  laser,  model  Powerlite  8010  by 
Continum  Corp.  ,  was  used  at  532  nm  and  800  mJ/pulse 
as  the  source  to  induce  Raman  effect.  The  repetition  rate 
for  this  laser  is  10  Hz  with  pulse  width  of  10  ns.  To 
reduce  possibility  of  the  plasma  formation  a  passive  pulse 
stretcher  was  designed  as  shown  in  Fig.  2.  The  laser 
beam  passes  through  three  adjustable  pinholes  (PHI, 
PH2,  and  PH3)  positioned  at  sufficient  distances  to  allow 
periodic  inspection  of  the  incident  beam  alignment. 
Because  the  laser  cavity  was  located  on  a  separate  optical 
table,  this  arrangement  was  found  to  be  a  very  effective 
method  for  realignment.  Also,  this  was  used  to  direct  a 
He-Ne  alignment  laser  beam  collinear  with  the  Nd:Yag 
for  initial  phase  of  the  alignment  procedure.  Wavelengths 
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Figure  2.  Schematic  diagram  of  the  laser  beam/sheet 
path  showing  the  delay  loops  and  the  alignment 
laser.  BS  (beam  |plitter),  PH  (pinhole),  M  (mirror), 
L  (spherical  lens),  and  CL  (cylindrical  lens). 


not  being  the  same,  this  initial  alignment  is  near  the  final 
one  which  can  ultimately  be  achieved  when  the  Nd:Yag 
laser  is  used.  This  was  found  to  be  a  safer  and  more 
efficient  approach  in  usage  of  the  Nd:Yag  laser  lifetime. 
In  Fig.  2,  the  laser  beam  faces  the  first  beam  splitter  at 
BS1  and  about  35%  of  its  intensity  is  split  towards  the 
M6  mirror.  The  remaining  goes  through  the  first  loop 
(B  S 1 -M4-B  S2-M5 -B  S 1 )  in  about  8  ns  of  time.  Note  that 
only  about  35%  of  this  remaining  energy  is  reflected  by 
the  mirror  M5,  the  rest  goes  through  the  loop  (BS2-M1 1- 
M10-M9-BS2)  in  another  8  ns  time.  Therefore,  the 
original  pulse  is  extended  in  time  into  three  major  pulses 
each  being  10  ns  in  duration  but  delayed  in  time  by  §hs 
and  16ns.  The  final  pulse  duration  is  about  36ns. 
Precision  laser  mirrors  with  damage  threshold  of  10  J/cm2 
(lO&s  pulse)  and  VIS  coating  from  Edmond  Scientific 
were  used.  The  laser  beam  is  shaped  by  passing  through  a 
Galilean  beam  expander  (LI  &  L2)  before  forming  a  thin 
sheet  via  two  cylindrical  lenses  CL1  and  CL2,  see  Fig.  2. 
The  entire  assembly,  consisting  of  mirror  M8  and  all  the 
lenses,  are  on  a  separate  optical  bench  which  can  be 
traversed  as  a  unit  vertically  (into  the  plane  of  Fig.  2)  as 
well  as  in  one  horizontal  direction.  This  gives  sufficient 
flexibility  to  adjust  and  fix  both  the  thickness  of  the 
sheet  and  its  vertical  and  horizontal  locations.  The 
vertical  and  horizontal  traverse  directions  provide  axial 
and  radial  traverse  capabilities,  respectively. 


The  scattered  Raman  signal  at  607  nm  was  detected  by 
a  nitrogen-cooled  256x256  CCD  from  Princeton 
Instrument  equipped  with  a  set  of  three  filters.  A  high- 
pass  OG570  filter  was  used  to  reduce  the  light  at  532hm 
by  less  than  10’5  and  at  607  by  not  more  than  0.9.  ,  A 
holographic  supemotch  filter  blocked  the  light  at  532nm 
better  than  10'4  while  less  than  0.8  at  607nm.  Finally,  a 
narrow  bandpass  interference  filter  at  610  nm  reduced 
the  scattered  light  at  532nm  to  better  than  10'4  while 
attenuating  by  about  0.5*  at  607.  Therefore,  an  optical 
density  of  better  than  13  was  achieved  at  53^nm  laser 
wavelength  with  the  Raman  signal  transmittance  of 
0.36.  This  arrangement  was  found  to  be  adequate  to 
ensure  removal  of  stray  light  while  passing  sufficient 
Raman  for  detection  at  607hm.  The  lens  package 
consisted  of  a  Nikon  PK-12/14  extension  ring  connected 
to  one  end  of  a  Noct-Nikkor  5^(nm  F/#  (1.2  to  16)  zoom 
lens  with  a  NIK-2734  close-up  lens  located  in  its  front. 

A  Summary  of  Previous  Results 

In  the  past  two  years,  results  from  injections  of  several 
fluids  into  an  ambient  under  both  sub-  and  supercritical 
conditions  at  sufficiently  high  Reynolds  numbers  to  be 
fully  turbulent  have  been  reported  in  the  same  test  facility 
used  here,  see  Chehroudi  et  al.  [8,9].  A  variety  of 
ambient  fluids  was  used  into  which  pure  N2  ,  He,  and  O2 
fluids  were  injected.  The  effects  of  chamber  pressure 
(density)  ranging  from  the  thermodynamic  subcritical  to 
supercritical  values  at  a  supercritical  chamber  temperature 
(based  on  critical  condition.  Pc,  Tc,  of  tire  injectant)  were 
observed  by  acquisition  of  shadow  images  from  the 
injector  exit  region  using  a  CCD  camera  illuminated  by 
short-duration  light  pulses.  At  sufficiently  low  subcritical 
chamber  pressures,  the  disturbances  on  the  jet  interface 
amplified  and  eventually  broke  up  downstream  into 
irregularly-shaped  small  entities.  Increasing  the  chamber 
pressure  caused  the  formation  of  many  small  ligaments 
and  droplets  at  the  interface  of  the  jet  only  within  a 
narrow  regime  below  the  thermodynamic  critical  pressure 
of  the  injected  pure  fluid,  resembling  a  second  wind- 
induced  liquid  jet  breakup.  At  even  higher  chamber 
pressures,  near  but  below  the  critical  pressure  of  the 
injectant,  the  expected  transition  into  a  full  atomization 
regime  to  produce  a  liquid  spray  was  inhibited  due  to 
sufficient  reduction  of  both  the  surface  tension  and  the 
heat  of  vaporization.  The  jet  appearance  changed  abruptly 
at  this  pressure  and  resembled  a  turbulent  gas  injection 
for  all  higher  chamber  pressures.  The  jet  initial  growth 
rate  was  plotted  together  with  available  data  on  liquid  fuel 
injection  in  diesel  engine  environments,  and  turbulent 
incompressible,  supersonic,  and  variable-density  jets  and 
mixing  layers.  The  resulting  plot  is  unique  on  its  own 
right  covering  four  order  of  magnitude  in  density  ratio, 
see  Fig.  3.  At  near-  and  super-critical  pressures,  these 
measurements  agreed  well  with  the  theoretical  growth 
rate  equations  proposed  by  Brown  [ 1 5]/Papamoschou 
and  Roshko  [16]  and  Dimotakis  [17]  for  incompressible 
but  variable-density  turbulent  mixing  layers.  This 
constitutes  the  first  quantitative  evidence  in  support  of 
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Jet  or  Spray  or  Mixing  Layer  Divergence  Angle  vs  Chamber-to-injectant  Density  Ratio 


Figure  3.  Initial  growth  rate  of  fluid  jets  as  a  function  of  the  chamber-to-injectant  density  ratio  for  variety  of 
cases.  iL. 


the  past  qualitative  observations  that  the  jet  appeared  to 
evolve  into  a  gas-like  behavior.  The  geometry  of  the  jet 
interface  had  also  been  examined  for  the  first  time  by 
fractal  analysis.  The  results  clearly  indicated  a  transition 
from  a  Euclidean  to  a  fractal  interface,  with  a  fractal 
dimension  close  to  values  measured  for  gaseous  turbulent 
jets.  This  provided  an  additional  quantitative  evidence  for 
the  hypothesis  that  the  jet  evolved  into  a  gas-like 
behavior.  Finally,  an  equation  was  proposed  based  on  a 
physical  model  proposing  that  at  the  point  of  transition 
from  liquid-like  to  gas-like  appearances  and  the  growth 
rates,  the  characteristic  time  of  the  vaporization  is  of  the 
same  order  as  the  interfacial  “bulge”  formation/separation 
time.  The  model  equation  agreed  well  with  the 
experimental  growth  rate  data. 

Experimental  Results  and  Discussion 

The  first  task  prior  to  any  image  acquisition  from  the  jets 
was  to  verify  basic  operation  of  the  entire  system.  For 
this,  it  was  decided  to  measure  the  Raman  intensity  from 
a  pure  gaseous  nitrogen  inside  the  high-pressure 
chamber  at  different  chamber  pressures.  The  ratios  of  the 
measured  Raman  intensities  (corrected  for  the  dark  count) 
should  be  proportional  to  the  density  ratios  calculated 
through  an  equation  of  state  or  thermodynamic  tables. 
Figure  4  shows  results  from  such  a  test. 


In  each  column  of  this  figure,  gray  scale  (or  false  color,  if 
seen  in  color)  representation  of  the  scattered  Raman  from 
the  laser  sheet  is  shown  along  with  its  5-row-averaged 
intensity  profile.  The  averaging  is  done  within  the  two 
thick  horizontal  lines  (covering  5  rows  of  pixels)  shown 
in  images.  The  lower  edge  of  the  laser  sheet  is  clearly 
seen  in  each  image  and  is  where  dark  black  zone  meets 
the  brighter  color  code.  After  removal  of  the  dark-count 
from  images,  the  measured  intensity  ratios  are  compared 
with  those  calculated  using  the  ideal  gas  law  for  nitrogen. 
The  agreement  was  quite  satisfactory  to  within  2  percent. 

When  no  LN2  was  injected,  the  image  acquisition 
procedure  for  each  set  consisted  of  an  initial  dark  count 
image  followed  by  images  at  three  different  chamber 
pressures  (6.9,  2.82,  and  1.46  MPa)  similar  to  Fig.  4. 
This  information  was  needed  to  remove  dark  count  level 
and  establish  reference  intensity  levels  at  the 
aforementioned  chamber  pressures.  They  also  served  as  a 
check  for  proper  operation  of  the  entire  system  setup. 
The  Raman  image  acquisition  from  the  LN2  jets  at  this 
same  chamber  pressures  started  subsequently. 

Typical  averaged  (40  frames)  raw  images  of  the  LN2  jet 
at  subcritical  and  supercritical  chamber  conditions  along 
with  images  characterizing  the  sheet  profile  at  the 
measurement  location  are  shown  in  Fig.  5.  The  actual 
measured  intensity  values  were  divided  by  2  to 
accommodate  the  plots  within  the  fixed  Y-axis  scale  of 
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6.90  MPa  (1000  psig) 


2.82  MPa  (400  psig) 


1.46MPa  (200psig) 


Figure  4.  Raman  intensity  measured  by  passing  the  laser  sheet  into  the  chamber  filled 
with  N2  at  three  different  pressures. 


the  software.  However,  the  conclusions  reached  from 
these  plots  are  unchanged  by  such  a  scaling.  Note  that  the 
Raman  intensity  profile  for  the  subcritical  case  is  higher 
in  the  center  of  the  jet  in  magnitude  and  appears  narrower 
in  shape.  With  LN2  injection,  it  was  observed  that  the 
measured  intensity  well  outside  the  jet  area  (representing 
the  chamber  density)  was  higher  than  when  injection  was 
turned  off.  This  is  due  to  the  interaction  of  the  laser  and 
the  jet,  being  particularly  severe  for  the  subcritical  case 
where  a  clear  sharp  interface  can  be  defined.  This 


interface  caused  reflection  of  a  fraction  of  the  laser 
contributing  to  higher  measured  Raman  intensity  values 
and  noise  levels.  As  an  example,  on  the  40-framed- 
averaged  images,  when  jet  flow  was  off,  the  pixel- 
averaged  measured  Raman  signals  within  a  fixed 
rectangular  region  where  laser  entered  the  chamber  were 
101.4  and  70.6  units  for  the  6.9  MPa  and  1.46  MPa  cases, 
respectively;  whereas  the  corresponding  values  for  when 
the  LN2  jet  was  flowing  were  116.4  and  98.0  units. 
Therefore,  this  defined  a  methodology  for  a  first  degree 


Supercritical  Subcritical  Laser  Sheet  Profile 

6.9  MPa  (1000  psig)  1.46  MPa  (200  psig)  6.9  MPa  (1000  psig) 


Figure  5.  typical  40-frame-averaged  false-color  images  of  the  jet  at  sub-and  supercritical  chamber  conditions. 
The  laser  sheet  profile  is  also  shown  to  be  reasonably  fat  within  the  measurement  region.  Measured  intensities 
are  shown  in  each  case. 


correction  to  the  raw  data  as  follows.  For  each  set  of  data 
at  a  given  condition,  Raman  images  were  collected  with 
no  injection  followed  by  the  Raman  images  from  the  jets 
as  described  earlier.  Using  the  example  above,  all  pixel 
values  for  the  6.9MPa  images  were  then  reduced  by 
(116.4  - 101.4  =  ih)  and  those  for  the  1.46  MPa  by  (98.0 
-  70.6=  27.4).  It  is  seen  that  this  laser/jet  interactional 


Figure  6.  40-frame-averaged  Raman  relative 
intensity  radial  profiles  at  two  different  distances 
from  the  injector  and  at  sub-  and  supercritical 
chamber  conditions. 

contribution  to  the  measured  raw  intensity  is  almost  twice 
as  much  for  the  subcritical  case  indicating  difficulty  in 
complete  removal  of  this  side  effect. 

The  test  conditions  for  the  Raman  imaging  were  selected 
to  correspond  and  complement  our  previously-reported 
shadowgraphy  studies.  Figure  6  Shows  typical  radial 
profiles  of  the  40-frame-averaged  Raman  intensities  as 
measured  with  reference  to  the  far  field  chamber  intensity 
values  at  two  distances  from  the  injector  exit  plan  and  at 
sub-  and  supercritical  chamber  pressures.  The  laser  sheet 
propagation  direction  is  from  right  to  left  in  this  figure. 
Limitation  for  the  number  of  frames  were  set  by  the 
combination  of  the  acquisition/processing  system  and 
the  practical  operating  time  of  the  facility.  However,  it 
was  found  to  be  reasonably  adequate.  It  is  seen  that  the 
centerline  intensity  surplus  (I-  Ich)  declines  away  from 
the  injector  as  well  as  when  chamber  pressure  is  raised  to 
a  supercritical  value.  Hence,  reasonable  trends  are 
observed.  One  point  of  concern,  however,  is  the  existence 
of  an  unusual  wide  tail  in  these  profiles  particularly  for 
the  subcritical  case.  This  issue  will  be  discussed  later. 
Figure  7  shows  similarity-type  plot  of  the  intensity 
surplus  at  different  axial  distances  and  at  sub-  and 
supercritical  conditions  along  with  a  proposed  model  for 


turbulent  jets.  Note  that  according  to  Wygnanski  and 


Figure  7.  Similarity-type  plot  for  the  measured 
Raman  intensity  surplus  (I  -  Ich)  at  sub-  (above) 
and  supercritical  (below)  chamber  conditions. 
R0.5  id  the  radius  at  which  (I  -  Ich)  drops  to  half 
of  its  maximum. 

Fiedler  [18]  fully  self-preservation  turbulent  jet  is 
observed  for  the  velocity  field  at  x/D  of  greater  than  40  in 
air  jet  when  Re  is  near  100,000.  So  et  al.  [19]  reported 
self-preservation  for  x/D  larger  than  about  20  for  binary 
gas  jet  at  Re  of  4,300.  Although  there  appear  to  be  some 
inconsistencies  in  this  criteria,  one  expects  to  reach  to  a 
similarity  profile  at  the  last  few  axial  distances  measured 
here.  What  emerges  as  immediately  obvious  in  Fig.  7  is 
the  wide  tail  in  the  profiles  particularly  for  the  subcritical 


Figure  8.  Similar  to  Fig.  7  but  results  by  other 
investigators  are  also  shown. 
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Figure  9.  Dark-count-corrected  Raman  intensity  (I) 
profiles  at  several  different  horizontal  traverse 
locations  but  at  a  fixed  distance  from  the  injector 
exit  (x/D=12.2).  The  Laser  sheet  is  parallel  to  the 
axis  of  the  jet  and  is  traversed  horizontally  towards 
or  away  from  the  camera.  The  symbol  “0”  means 
laser  sheet  passes  through  the  jet  axis;  “lccw” 
means  a  0.5mm  shift  from  the  jet  axis  towards  the 
camera  direction;  and  “lew”  means  a  0.5mm  shift 
from  the  jet  axis  away  from  the  camera.  Pr= 
0.43(top),  0.8 3 (middle),  and  2.03(bottom). 


case.  Also,  the  profiles  appear  to  be  flatter  near  the 
center  for  the  subcritical  case  exhibiting  a  tendency  for  a 


r/R0.5 

Figure  10.  Normalized  plot  of  the  Modified  and 
center  (i.e.  “0”)  profiles  shown  in  Fig.  8.  x/D=12.2 


near  top-hat  shape  in  this  region.  In  Fig.  8  results  by  So 
et  al  [19]  (gaseous  jets)  and  Oschwald  et  al  [13] 
(supercritical  cryogenic  jets)  are  also  shown.  They  clearly 
point  to  difficulties  in  obtaining  high  spatial  resolution 
and  fidelity  information  at  these  distances  from  the 
injector  exit  plane.  Note  that  the  outer  part  of  the  So  et  al. 
profiles  agree  with  our  data  whereas  those  by  Oschwald  et 
al.  follows  the  model  at  x/D  of  8.4.  Although  only  two 
profiles  are  shown,  data  by  Oschwald  et  al.  [13]  appears 
to  be  more  intuitive  as  it  suggests  a  gradual  change  of  the 
shape  with  distance  towards  the  far  field  model. 

To  investigate  the  causes  for  the  aforementioned  profile 
broadening,  the  vertical  laser  sheet  was  traversed 
horizontally  towards  and  away  from  the  camera  and 
results  were  analyzed  at  several  selected  axial  positions. 
Figure  9  shows  results  from  such  an  experiment.  A 
modified  profile  is  then  constructed  using  the  maximum 
values  of  the  Raman  intensity  profiles  at  each 
horizontally-traverse  position  as  indicated  in  Fig.  9.  It  is 
seen  that  this  modified  profile  is  narrower  than  the  profile 
taken  when  the  laser  sheet  slices  the  jet  containing  the 
centerline  (i.erjthe  “0”  case  in  Fig.  9).  This  effect  is 
particularly  ^strong  at  the  subcritical  conditions  and 
almost  negligible  when  supercritical.  This  suggests  that 
reflection  from  the  jet  boundary  may  induce  additional 
Raman  effect  in  this  zone  causing  broadening  of  die 
intensity  surplus  profiles  indicated  earlier.  Considering 
the  injector  hole  diameter  used  in  this  study,  such  a 
broadening  is  to  be  expected. 

Figure  10  shows  similarity-type  plots  for  the  modified  and 
their  corresponding  center  (identified  as  “0”  in  Fig.  9) 
profiles.  The  R0.5  point  is  determined  for  each  profile 
independently.  The  single  solid  circle  shows  the  end  of 
the  1%  point  to  be  at  about  2.5  times  the  full  width  at  the 
half  maximum  (FWHM)  value  (i.e^-2*  R0.5).  These 
modifies  plots  suggest  lack  of  type  similarity  seen  in 
gaseous  turbulent  jets  at  the  lowest  subcritical  condition 
(i.^:Pr=0.43).  In  contrast,  tendency  towards  such  a 
siitmarity  is  observed  at  Pr  of  0.83  and  2.03.  Admittedly, 
the  agreement  is  not  perfect.  Comparing  the  modified 
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instead  the  FWHM  appears  to  be  a 
safer  parameter  to  determine.  Also, 
there  is  doubt  if  density  values  can  be 
measured  to  within  1%  accuracy. 
Nevertheless,  Brown  and  Roshko  [20] 
reported  that  their  1%  thickness 
definition  corresponded  fairly  well  with 
the  extent  of  the  mixing  visible  on  the 
shadowgraph  images.  However,  the 
spread  angle  based  on  the  density 
profile  was  always  larger  than  that  on 
the  velocity  profile. 

The  FWHM  at  each  axial  position  was 
Figure  1 1.  Normalized  FWHM  of  the  density  surplus  radial  profiles  as  a  determined  and  results  along  with  data 
function  of  the  normalized  Distance  from  the  injector.  by  other  investigators  are  shown  in  Fig. 


x/b:  Re: 

Osehwald  et  al 

K-  MPa-  K; . : . Ratio- . ram: . . 

N2/N2 . V  Vis:  4:  298;  1.i?j . 3.34; . 

8.42:  ilSOOO  in  Figs.  7.  8 

Oschwald  et  al 

N2/N2  140;  4:  298;  J.17: . 12 'A' . 

Chehroudi  et  ai 
Chehroudi  et  ai 

N2/N2 . 1 . 95 : . 6.9 : "  ‘  2*95: . j.A . . .  ?.??. 

N2/N2 . . l!st*"295: . 0*4:* . 40.6;  0.5051  100 

4.8  to  24.4  3500*6: in  Figs.  8,1*6, ii 

4.8  to  24.5.  12000  in  Figs.  10.11 

So  et.  ai. 

. — •!■ . 

5.1:  5 666: in  Figs.  8,  11 

6.4:  5000 : in  Figs.  8.  ii 

So  et.  al. 

Richards  &  Pitts 
Richanis  &  Pitts 

26-80  :  4000  in  Fig.  11 

40-1 20 :  25666: in  Fig.  11 

Table  L  Some  information  extracted  from  works  by 
other  investigators  reported  here. 


Qiamber/lnjectant  Density  Ratio 


11.  More  information  on  the 
conditions  under  which 
experiments  were  performed  is 
given  in  Table  1.  Note  that  except 
for  Chehroudi  et  al.  [8,  9]  and 
Oschwaldetal.  [13]  all  others  are 
for  injection-  of  gaseous  fluids 
into  an  ambient  gas  at  subcritical 
pressure  (based  on  the  injectant 
critical  pressure).  The  results  of 
present  work  in  Fig.  11  show  an 
increasingly  larger  spreading  rate  as  chamber  pressure  is 
raised.  The  data  at  the  supercritical  condition  (Pr  =  2.03) 
approaches  that  by  Richards  and  Pitts  [21]  acquired  at 
density  ratio  of  larger  than  1  even  though  it  is 
substantially  smaller  than  ours.  Unfortunately,  for  our 
jets,  data  at  larger  distances  are  not  available  at  present 
time  to  develop  a  more  comprehensive  comparison 
between  the  cases.  Richards  and  Pitts  [21]  conclude  that 
if  care  is  exercised  to  ensure  that  flow  is  free  of  buoyancy 
and  coflow  effects,  the  spreading  rate  in  variable-density 
jets  is  independent  of  the  initial  density  ratio,  velocity 
profile,  turbulence  level,  and  conforms  with  the  results  of 
constant-density  by  others  and,  in  addition,  they 
propose  a  slope  in  the  range  between  0.212  to  0.220  for 
this  linear  jet  growth  rate  equation,  see  Fig.  11.  Our  data 
at  Pr=2.03  gives  a  least-square  fit  to  a  line  with  a  slope  of 
0. 102,  almost  half  of  that  by  Richards  and  Pitts  [  13],  It  is 
therefore  possible  that  at  some  high  enough  injectant-to- 
chamber  density  ratio  the  spread  rate  universality 
indicated  by  them  breaks  down  and  one  observes  a 
retarded  growth  rate  for  variable-density  turbulent  jets. 


Figure  12.  Comparison  of  the  spreading  angle  measured 
using  shadowgraph  and  Raman  techniques  using  twice 
the  FWHM  values.^ 

profiles  with  the  ones  taken  at  the  center,  one  can  deduce 
the  FWHM  from  these  center  profiles  with  not  more 
than  a  5%  inaccuracy  despite  the  existence  of  the 
broadening  effect.  When  available,  the  best  value  for  the 
FWHM  were  considered.  The  observed  broadening  in 
this  work  prevents  us  from  using  the  1%  free  stream 
density  definition  used  by  Brown  and  Roshko  [20]  and 


The  measured  spreading  rate  using  Raman  is  compared 
with  the  spreading  rate  we  reported  in  the  past  from  the 
shadowgraph  visualizations.  The  results  are  plotted  in  Fig. 
12.  It  can  be  seen  that  the  spreading  rate  at  twice  the 
FWHM  values  agrees  well  with  the  shadowgraph 
measurements.  At  subcritical,  a  factor  less  than  2  would 
give  a  better  agreement  as  implied  from  Figs.  10  and  12. 
Thus  the  correspondence  of  the  two  different  techniques 
for  measuring  the  width  of  the  jet  has  been  established. 
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Conclusions  and  Summary 

Two-dimensional  Raman  signal  scattered  from  a 
cryogenic  LN2  jet  was  acquired  at  several  distances  from 
the  injector  both  at  sub-  and  supercritical  conditions. 
Reflection  of  the  laser  off  of  the  jet  boundary  produced 
additional  Raman  effect  in  this  zone  causing  profile 
broadening.  This  effect  was  demonstrated  through 
systematic  horizontal  traverses  of  the  laser  sheet  towards 
and  away  from  the  camera.  Within  experimental 
accuracy,  the  modified  intensity  surplus  plots  indicated  a 
similarity  behavior  of  the  type  observed  for  classical 
constant  and  variable-density  gaseous/liquid  jets  only  at 
near-  and  super-critical  conditions.  Because  of  this 
broadening  effect  and  as  far  as  a  good  measure  of  the  jet 
thickness  and  growth  rate  are  concerned,  the  FWHM  of 
the  profile  is  a  better  estimator  to  use  than  the  1%  (above 
the)  far  field  value.  The  growth  rate  of  the  LN2  jet,  using 
slope  of  the  linear  FWHM/D  versus  x/D  plot,  increases 
monotonically  as  chamber  pressure  reaches  critical  and 
supercritical  values.  A  tendency  exits  for  the  above 
linear  FWHM  plot  at  supercritical  condition  to  approach 
those  reported  for  gaseous  jets  affirming  our  gas-like 
appearance  reported  in  the  past.  It  is  also  shown  that  a 
plot  of  twice  the  FWHM  versus  density  ratio  agrees 
with  the  jet  thickness  measured  through  our 
shadowgraphy  studies.  This  establishes  an  equivalency 
between  these  two  different  jet  thickness  determination 
methods. 
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